Effect of the tensor force on the charge-exchange spin-dipole excitations of ^^Pb 



O 






> 

in 

o 
o 



% 



ex. Baii'2), H.Q. Zhangi^2)^ jj. Sagawa^), X.Z. Zhang^), G. Colb'^\ and F.R. Xu^' 

^' China Institute of Atomic Energy, China and 

^' School of Physics and State Key Laboratory of Nuclear Physics and Technology, Peking University, China 

'^' Center for Mathematics and Physics, University of Aizu, Aizu-Wakamatsu, Fukushima 965-8560, Japan 

'^' Dipartimento di Fisica, Universith degli Studi di Milano and IN FN, Sezione di Milano, 20133 Milano, Italy 

The charge-exchange spin-dipole (SD) excitations of ^"^Pb are studied by using a fully self- 
consistent Skyrme Hartree-Fock plus Random Phase Approximation (HF-I-RPA) formalism which 
includes the tensor interaction. It is found, for the first time, that the tensor correlations have 
a unique, multipole-dependent effect on the SD excitations, that is, they produce softening of 1~ 
states, but hardening of 0" and 2~ states. This paves the way to a clear assessment of the strength 
of the tensor terms. We compare our results with a recent measurement, showing that our choice of 
tensor terms improves the agreement with experiment. The robustness of our results is supported 
by the analytic form of the tensor matrix elements. 
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The nuclear effective interactions like the zero-range 
Skyrme forces have been quite successful to describe 
many nuclear properties. These forces are fitted using 
empirical properties of uniform nuclear matter, together 
with masses and charge radii of selected reference nu- 
clei. They describe in a reasonable vifay the global trends 
of the ground-state properties along the nuclear chart 
(binding energies, radii and deformations). Properties of 
excited states such as vibrations and rotations have been 
studied successfully as well, allowing a large amount of 
physical insight [l|, l3| • In the quest for a universal local 
Energy Density Functional (EDF) for nuclei, the Skyrme 
framework is often used as a starting point. 

While most of the Skyrme parameter sets which have 
been widely used are purely central, many groups have 
recently devoted attention to the role played by the zero- 
range tensor terms that can be added (see Refs. [j- 
Il0l|l. This blooming of theoretical studies has followed 
the claim by the authors of Ref. [ll|, that the tensor 
force is crucial for the understanding of the evolution of 
the single-particle energies in exotic nuclei. 

There exist different strategies to fix the tensor part of 
the interaction. One can be inspired by a bare or a G- 
matrix interaction 0, [l^l ■ Since the tensor force affects 
the spin-orbit splittings as described below, another pos- 
sibility is to add it to an existing Skyrme set and try to 
reproduce at best the evolution of single-particle states 
along isotopic or isotonic chains [a, |g|- At present, the 
most accurate and systematic way to produce effective 
interactions with the tensor, has been a full variational 
procedure to fit the tensor and the central terms on equal 
footing [31. All these attempts have produced results 
which are not at all conclusive. 

The studies based on a refit of all Skyrme parame- 
ters suffer from the drawback that the tensor force has 
a moderate effect on the ground state quantities such 
as the total binding energies. The single-particle ener- 
gies do not lie within the EDF framework and can be 
affected by correlations like particle-vibration coupling. 
Then, we follow in this work the idea that collective ex- 



citations (especially the spin-dependent ones) may be a 
better candidate to constrain the tensor force. 

Recently, self-consistent HF-(-RPA schemes with tensor 
interactions have been developed |13l - [l5| . The Gamow- 
Teller (GT) and charge-exchange 1"*" spin-quadrupole 
(SQ) transitions in ^"Zr and ^°^Pb have been studied 
in Refs. [Ij, [lj|, whereas the non charge-exchange mul- 
tipole responses of several magic nuclei have been calcu- 
lated in Ref. [15l|. In this letter, we study the charge- 
exchange spin-dipole (SD) excitations of ^°*Pb, inspired 
by recent accurate measurements [lg|. Our specific goal 
is to find a clear, unambiguous signature of the effec- 
tive tensor force. Spin-isospin collective modes have been 
instrumental for the understanding of nuclear structure 
since almost three decades [ij, ]1^. In particular the 
total SD strength distribution in ^°Zr has been mea- 
sured [i^, [20] and this has allowed to extract conclu- 
sions on the neutron radius |21| . Neutron radii have been 
shown to be strongly correlated with the features of the 
neutron matter equation of state (EOS) which, in turn, 
has relevant implications for neutron stars [22-24,]. How- 
ever, for the determination of the effect of the tensor force 
which is strongly spin-dependent one needs to know sep- 
arately the strength distributions of the J'^ = 0^, 1^ and 
2^ components. 

The triplet-even and triplet-odd zero- rang e tensor 
terms of the Skyrme force are expressed as [l^, [2^ 
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In the above expression, the operator k = (Vi — V2) /2i 
acts on the right and k' = — (V^ — V2) /2i acts on 
the left. The coupling constants T and U denote the 



strengths of the triplet-even and triplet-odd tensor inter- 
actions, respectively. 

The main effect of the tensor terms on the HF results 
is a modification of the spin-orbit potential, which reads 



The model independent sum rule which is known to hold 
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In this expression, g=0(l) labels neutrons (protons), pq 
are the densities while Jq are the so-called spin-orbit den- 
sities. Their definitions can be found in Refs. |3l-ll0l [l3| . 
The first term in the r.h.s of Eq. 1^ comes from the 
Skyrme two-body spin-orbit interaction, whereas the sec- 
ond term includes both a central exchange and a tensor 
contribution, that is, a — ac+ctr and /? = (5c+Pt- their 
complete expressions can be found e.g. in Ref. [3,[3]. In 
this letter, we employ different Skyrme parameter sets 
which include the tensor terms. The sets T/J have been 
introduced in Ref. [7|. The set SLy5-|-Tu, is a set in 
which the tensor terms are added perturbatively on top 
of the existing force SLy5: the tensor parameters T and 
U are obtained by the low-g limit the G-matrix calcu- 
lations ^2j]. One should notice that the tensor part of 
SLy5-|-T.,„ is different from that of SLy5-|-T which was 
adopted in Ref. [5]. We tried to use the set SLy5-|-T, but 
its effect on SD states is larger and do not match the SD 
experimental data in ^°*Pb [16|. The values of T, U, a 
and /? for the adopted interactions are listed in Table HI 



TABLE I: The tensor strength parameters T and U of Eq. 
((l| as well as the a and /? values of Eq. ((2)1 . All values are in 
MeV-fm^ 



Force 


T 


U 


CtT 


Pt 


etc 


Pc 


SLy5+T„ 


820.0 


323.4 


134.76 


238.2 


80.2 


-48.9 


TU 


259.0 


-342.8 


-142.8 


-17.5 


82.8 


-42.5 


T22 


356.0 


-217.5 


-90.6 


28.9 


90.6 


-28.9 


T32 


613.1 


-2.3 


-96.5 


79.5 


96.5 


-19.5 


T43 


590.6 


-147.5 


-61.5 


92.3 


121.5 


27.7 


T44 


520.1 


21.5 


9.0 


112.8 


111.0 


7.1 



In our model, we first solve the HF equations in co- 
ordinate space. The unoccupied levels are found by di- 
agonalizing the mean field in a harmonic oscillator basis 
(up to a maximum value of the major quantum number 
Nmax — 12). We then perform fully self-consistent RPA 
by including both the two-body spin-orbit and tensor in- 
teractions. It has been checked that the adopted basis is 
large enough to make the results stable. 

The charge- exchange SD operator is defined as 



oi = E4^'W^^)^ 



(3) 



The n— th energy weighted sum rules m„ for the A-pole 
SD operator are defined as 



:(i±) = E^ri(*io±io)i^ 



(4) 



ism^(i_)-m^(<+) 
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Z(r2)p). This sum 



rule has been shown to be fulfilled with 1% accuracy by 
the numerical calculations. 



TABLE II: The SD sum rules mo and mi for ^"**Pb with 
and without the tensor terms. A_E is the difference between 
mi /mo calculated with and without tensor. 







without tensor 


with tensor 




force 


A" 


mo 


m,i 


mi/mo 


mo 


mi 


mi/mo 


AS 




0" 


158.6 


4718 


29.7 


171.9 


5168 


29.9 


0.2 


SLy5 


1" 


432.0 


11746 


27.2 


440.0 


10111 


23.0 


-4.2 


+T„ 


2- 


646.0 


13742 


21.3 


657.4 


14408 


21.3 







sum 


1236.4 


30206 


24.4 


1269. 


29256 


23.0 


-1.4 




0- 


158.6 


4559 


28.8 


163.4 


6074 


37.2 


8.4 


TU 


1" 


434.5 


11581 


26.7 


431.0 


10376 


24.1 


-2.6 




2- 


657.3 


13710 


21.2 


646.2 


14296 


22.1 


1.0 




sum 


1240. 


29850 


24.1 


1241. 


30746 


24.8 


0.7 




0" 


157.6 


4650 


29.5 


163.7 


5943 


36.3 


6.8 


T22 


1" 


434.5 


11771 


27.1 


433.9 


10504 


24.2 


-2.9 




2" 


645.9 


13812 


21.4 


647.1 


14327 


22.1 


0.7 




sum 


1238.0 


30233 


24.4 


1244.7 


30774 


24.7 


0.3 




0" 


157.2 


4698 


29.9 


166.9 


6479 


38.8 


8.9 


T32 


1" 


435.3 


11886 


27.3 


435.3 


9928 


22.7 


-4.6 




2' 


646.1 


13897 


21.5 


649.6 


14619 


22.5 


1.0 




sum 


1238.6 


30481 


24.6 


1253.2 


31026 


24.8 


0.2 




0" 


154.8 


4693 


30.3 


164.0 


6170 


37.5 


7.2 


T43 


1" 


440.3 


12138 


27.6 


444.1 


10366 


23.3 


-4.3 




2" 


645.5 


14067 


21.8 


649.4 


14675 


22.6 


0.8 




sum 


1241. 


30898 


24.9 


1258. 


31211 


24.8 


-0.1 




0~ 


155.6 


4811 


30.9 


163.2 


5637 


34.5 


3.6 


T44 


1" 


436.3 


12174 


27.9 


440.6 


10854 


24.6 


-3.3 




2- 


664.2 


14059 


21.8 


649 


14406 


22.2 


0.4 




sum 


1256. 


31044 


24.7 


1253. 


30897 


24.7 


0.0 



We performed two kinds of calculations. In the first 
one, the tensor terms are neither included in HF nor 
in RPA. In the second one, the tensor terms are in- 
cluded both in HF and in RPA. The terms containing 
the spin-orbit densities Jq which arise from the central 
momentum-dependent part of the Skyrme interaction are 
included in all the following HF and RPA calculations. 
Previously, it had been found that the effect of tensor 
correlations in HF is large for the Gamow- Teller mode 
and also, to some extent, for some low-lying non charge- 
exchange excitations like the first 2+. This is largely due 
to the fact that unperturbed GT transitions are exactly 
those among spin-orbit partners. The unperturbed p-h 
energy of the lowest 2+ is also determined by the spin- 
orbit splitting. On the other hand, this is not the case 
for the SD excitations: the average unperturbed energies 
are not much affected by the spin-orbit splittings since 
they are the Ifvjj type excitations. 

The numerical results of the HF+RPA calculations 
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FIG. 1: (Color online) Charge-exchange SD_ strength distri- 
butions in ^''^Pb. In the panels (a), (b), and (c) the RPA re- 
sults obtained by employing the interaction SLy54-Tuj for the 
multipoles 0~ , 1~, 2~ are displayed. In panel (d) we show 
the total strength distribution. Panels (e), (f), (g) and (h) 
correspond to similar results when the parameter set T43 is 
employed. All these discrete RPA results have been smoothed 
by using a Lorentzian averaging with a width of 2 MeV and 
compared with experimental findings. The excitation energy 
is with respect to the ground state of '^"^Bi. The experimen- 
tal data are taken from Ref. ]W^ - See the text for details and 
discussion. 



shown. From Fig. [TJe), we see that the calculated 0~ 
strength is concentrated in one peak which is shifted up- 
wards by about 1.3 MeV by the tensor correlations. In 
FiglUf), the RPA tensor correlations move the 1" peak 
downwards and split it into three peaks, in qualitative 
agreement with the bump-like experimental strength. In 
the case of the 2" [Figdjg)], the main peaks in the high 
energy region are rather near to the experimental main 
peak. Therefore, as shown in Figlljh), the inclusion of 
the tensor terms in HF-I-RPA can make the calculated 
main peak of the total SD strength coincide with the 
main measured peak. However, in low energy region 
the agreement is not good compared with the experi- 
mental result in the case of SLy5-|-Tu, (due also to the 
fact that with tensor a peak appears below the g.s. of 
the daughter nucleus). The spin-orbit density gives con- 
tributions to the p-h matrix elements of spin-dependent 
excitations [26]. We find out in the case of T22 param- 
eter set that the two-body interaction generated by the 
terms containing Jq from the central exchange interac- 
tion have a much smaller effect on the SD excitations 
than the two-body tensor interaction |27| . 

We would like at this stage to obtain a better under- 
standing of this peculiar role of tensor interactions. The 
diagonal matrix element of the triplet-even (TE) term on 
a state with multipolarity A can be expressed as 111 



V, 



w 



TE 




in terms of the reduced matrix elements of the operator 
Ofe.A = E^W^ «> (V^ «) y^(z))«](^ and 6j symbols. In 
Eq. (O, the notation k = \/2k + 1 is used. For the SD 
excitations, taking £ = and k=k'=l, Eq. ([5|) gives 



with the forces T43 and SLy5-|-Ttu are shown in Fig. 
[TJ They are compared with experimental data obtained 
by Distorted Wave Impulse Approximation (DWIA) and 
multipole decomposition analysis of the (p,n) reaction 
|la |. From Figdja) and (b) one can see that in the case 
of the T43 interaction the main peaks of the 0~ and 1~ 
strength distributions are shifted upwards by about 7.5 
MeV and downwards by about 5 MeV, respectively, due 
to the tensor correlations. There are several 2~ peaks (cf. 
FigHl^c)). The peak at excitation energy E^ «17.7 MeV 
is moved upward by about 2 MeV by including tensor, 
and comes close to an experimental peak, while another 
peak at E^ «3.9 MeV is shifted downwards by about 0.6 
MeV and is also eventually close to the observed low en- 
ergy peak. For the total SD strength in Fig. [TJd), it is 
remarkable that the main peak at 26 MeV is shifted to 
21 MeV when tensor is included, and this provides good 
agreement with the experimental data. 

In the same figure, the SD strength distributions in 
^°^Pb calculated by using the set SLy5+Tt„ are also 



1 



0- 



V, 



(A) 



TE 



T7;T< -1/6 } \{p\\d^,x\\h)\' for A 



12 



1/50 



(6) 
The TO tensor part is also expressed in a similar way as 



^TO - ^2 




|(p||di,A||/i)P for A = 




We can see in Eqs. ([6]) and (O that the diagonal p-h 
matrix element in the 0~ case is the largest, and that for 
1~ is the next. The effect on 2~ is rather small. It should 
be noticed that these relative strengths of the Skyrme 
tensor interactions on each multipole are similar to those 
obtained from the finite-range tensor interactions both in 
magnitude and in sign [28[ . We can sum the TE and TO 
direct matrix elements as 



V, 



w 



Vr, 



(A) 



TE 



Vn 



(A) 



TO 



axT + hxU. 



(8) 



The proper antisymmetrization is easy to obtain for con- 
tact interactions and gives, in the isovector channel, 



V, 



w 



1 



1. 



T,AS 



= h77«AT+-foAC/](ri-r2). 



(9) 



Since the coupHng constant T is positive for the interac- 
tions we considered, Vj.^ is repulsive for the 0^ and 2^ 
case, while it is attractive for 1^ . The Vj,q part may 

contribute with the same sign as the Vj,J: one if the value 
of U is negative. For the T/J family, the value of U is 
negative or small positive, so that the Vj,q contributions 
have the same multipole dependence or almost negligible. 
All together, the tensor correlations are strongly repul- 
sive for 0~ and weakly repulsive for 2~ in general. For 
1~, the net effect will be attractive. For SLy5-|-Tiu, the 
value of U is positive and will give opposite contributions 
to those of T. However, the T value is much larger than 
the value of U so that the same argument given for the 
TJJ family will hold. One can see from Table HIl that Eq. 
^ provides a very effective guideline for interpreting the 
numerical results of microscopic RPA. In the same Table 
we also provide values of the sum rules mo and mi for 
the different multipoles and effective forces. 

In the parameter sets T43 and T44 used for Table I, 
the spin-orbit strengths are larger than in the other pa- 
rameter sets. However, these large spin-orbit strengths 
cancel substantially with the effect of the terms a and /3 
of Eq. (2), and the net results are quite similar for the 
single particle spectra of the parameter sets in Table I. 

In summary, we have studied the effect of tensor in- 
teractions on the charge-exchange t_ SD excitations of 
208pj^ by using the self-consistent HF-I-RPA model. We 
have demonstrated clearly for the first time that tensor 



correlations have a specific multipole dependence, that 
is, they produce a strong hardening effect on the O~mode 
and a softening effect on the 1~ mode. A weak hardening 
effect is also observed on the 2~ mode. The character- 
istic effect of the tensor force was further examined by 
using the analytic formulas based on the multipole ex- 
pansion of the contact tensor interaction. It is found 
that both the TE and TO tensor interactions do provide 
the characteristic multipole dependent effect on the SD 
excitations. These effects are shown to be robust irre- 
spective of the adopted Skyrme force with tensor inter- 
actions. Our calculated results are compared with recent 
SD excitation spectra obtained in the (p, n) reaction on 
the ^'^^Pb target. The softening and the hardening effects 
produced by tensor correlations on the 1~ and 2~ modes 
are confirmed by comparing the experimental data and 
the calculated results with and without tensor correla- 
tions. Consequently, from the study of SD excitations, 
we are able to give a clear constraint on the effective 
tensor force which was missing so far. 
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